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ABSTRACT 

We present observations of M31LRN 2015 (MASTER OT J004207.99+40550I.I), discovered in M31 
in 2015 January, and identified as a rare and enigmatic luminous red nova (LRN). Spectroscopic and 
photometric observations obtained by the Liverpool Telescope showed the LRN becoming extremely 
red as it faded from its My = —9.4±0.2 peak. Early spectra showed strong Ha emission that weakened 
over time as a number of absorption features appeared, including Na i D and Ba ii. At later times 
strong TiO absorption bands were also seen. A search of archival Hubble Space Telescope data revealed 
a luminous red source to be the likely progenitor system, with pre-outburst Ha emission also detected 
in ground-based data. The outburst of M31LRN 2015 shows many similarities, both spectroscopically 
and photometrically, with that of V838 Mon, the best studied LRN. We finally discuss the possible 
progenitor scenarios. 

Keywords: binaries: general — novae, cataclysmic variables — stars: individual (M31LRN 2015, V838 
Mon) — stars: peculiar — supergiants 


1. INTRODUCTION 

Luminous red novae (LRNe) are a rare type of stellar 
transient typically more luminous than most novae, but 
fainter than supernovae (SNe), and characterized by a 
reddening color as they fade. The canonical system is 
V838 Monocerotis, which produce d an outburst in early 
2002, peaking at My = —9.95 (Munari et al. 2002b|) 
and l ater generating a spectacular light echo (|Bond et al.' 
20031. The only extragalactic r ed transient with a simi¬ 
lar peak magnitude is M31-RV (Rich et al.||1989|). Addi¬ 
tionally, several other extragalactic transients show some 
similar properties to LRNe, but are sig nificantly more 
lumin ous at pea k (e.g., M85 OT200 6-1, Kulkarni et al. 
2007| PTE 10fqs,|Kasliwal et al.|2011|), and often grouped 
as intermediate luminosity red transients (ILRTs). 

A number of mechanisms have been proposed to ex¬ 
plain these events, inclu ding (1) a classical nova mecha¬ 
nism ([Shara et al.|| 2010|, (2 ) an unusual SN mechanism 
(e.g. i Lovegrove fe Woosley||2 013|), (3) giant plane t cap¬ 
ture dRetter &: Marom|j2005 ), (4) stellar mergers ( Soker 


nova in M31. 


fc Tyfenda]|20t 

son et al. 1120091 


200j and 

A 


5) extreme AGB stars (Thomp- 


Kochanek (2011), in a study ot ILKi's 
yiN 2008b and NGG 300 (JT'2008, suggested that LRNe 
and ILRTs may have different physical origins, favoring 
the ext reme AGB scenario a s the likely mechanism for 
ILRTs. Tylenda et al. (2011) presented compelling evi¬ 
dence that V1309 Scorpii, a Galactic LRN, was a merger 
of a contact binary. 

MASTER OT J004207.99-k405501.1 (hereafter 
M31LRN 2015) was discovered on 2015 January 
13.63 UT by the MASTER-Kislovodsk auto-detection 
system at 00*'42“07".99 -f40°55'01 "l (J2000) with 
an unfiltered magnitude of 19.0 (Shumkov et al. 


2015b). The transient then brightened significantly from 


January 14.71 ( 

Ovcharov et al. 

2015) 

taken by 

Uvcharov et al. 


on 2 


2015). A spectrum was 


which revealed the presence ot Ha emission, leading to 
the tentative classification of the transient as a classical 


have been reported bv 

Dong et al. 

(2015) 

, Wagner et al. 

(2015 

1, and Adams et al. (|2U15|). 



As part ot our extensive optical tollow-up campaign 

for M31 novae (e.g. Shatter et al. 

2011 

Williams et al. 


spectroscopically with the Liverpool Telescope (LT), a 
fully robotic 2m telescope on the island of La Palma. 

2. PHOTOMETRIC OBSERVATIONS 

We obtained multi-color photometric observations of 
M31LRN 2015 usi ng the 10:0 CCD camera on the LT 


(Steele et al. 
taken in eac 


2004). Typically 3 x 120 s exposures were 
filter. 


which were decreased to 3 x 60 s 
at later times as M31 became more difficult to ob¬ 
serve. The B, V and i'-band photometric data we re cal¬ 
ibrated using eight stars from Massey et al. (2006), with 
the z' -band magni t udes c omputed using transtorniations 
from Jordi et al. (2006). The z'-band data were cali- 


brated using a single SDSS DR9 (Ahn et al. 2012) star 
(the only star in the field suitable for the calibration), 
J004225.86-1-405651.4, which may have affected the reli¬ 
ability, but by comparison to the i'-band photometry, we 
determined any systematic effect is small (< 0.1 mag). 
The photometry is presented in Table and a light curve 
combining published data with those taken with the LT 
is shown in Figure 

The most striking feature of the fight curve is the in¬ 
creasingly red color after peak, from {B — i') = 0.650 ± 
0.007 on January 19.89 to {B — i') = 3.54 ± 0.09 on 
March 1.83. Our data imply optical peak was reached be¬ 
tween 2015 January 19.89 and January 22.84. M31 has 
a distance modulus of (m—iH)o = 24.43 ± 0.06 (|Freed- 
man & Madore 1990). If we assume foreground red¬ 
dening ol = 0.062 (Schlegel et al. 1998) and 

maximum total extinction at that p osition in M31 of 
E(^B-y) = 0.18 (Montalto et al.||2009), we conservatively 
assume that M31LRN 2015 was subject to reddening of 
Ei^b-v) = 0.12±0.06. As such, the absolute peak magni¬ 
tude of the transient was My = —9.4 ±0.2. It is possible 
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Table 1 

LT photometry of M31LRN 2015. 


Date (UT) 


2015 Jan 17.87 
2015 Jan 19.89 
2015 Jan 22.85 
2015 Jan 27.85 
2015 Jan 30.88 
2015 Feb 14.83 
2015 Feb 26.84 
2015 Mar 1.83 
2015 Mar 6.83 
2015 Mar 9.84 
2015 Mar 12.83 


B 


16.172 ±0.003 
15.954 ±0.004 
15.965 ±0.003 
17.013 ±0.008 
17.431 ± 0.009 
18.478 ±0.029 
19.178 ±0.041 
19.243 ± 0.088 


V 


15.657 ± 0.002 
15.459 ± 0.003 
15.403 ± 0.002 
16.080 ± 0.005 
16.303 ± 0.004 
16.720 ± 0.004 
17.260 ±0.011 

17.782 ± 0.036 


15.463 ±0.006 
15.304 ±0.006 
15.198 ±0.005 
15.495 ±0.006 
15.542 ±0.006 
15.485 ±0.006 
15.649 ± 0.006 
15.704 ±0.009 
15.854 ±0.014 
16.004 ±0.017 
16.121 ± 0.025 


15.230 ±0.009 
15.276 ±0.015 
15.248 ±0.015 
15.361 ± 0.023 


Note. — BVi' calibrated against stars from |Massey et al.| (|2006f, z' calibrated 
against star from SDSS. 



Figure 1. Light curve of M31LRN 2 015 from LT obse r vation s 
ffilled circles), as well as data from lOvcharov ah] ||2015||. 
Kurten^gyet^l.l | |2015a|bl l, [Shumkov et al.] | |2U15aj | and [l‘abrika| 


et al.|(|20i5| (all open circlesj. it-band data is represented by blue 

r\riirt1-o- 1/ rrvocirf T/T/ riT’ano'ci- F? raA- T lifrl-i+ fTTcnr- Aarlr frroTr- 


points; V, green; W, orange; R, red; /, light gray; i', dark gray; 
z', black. The light blue and green lines represent the B and V 
photo r netry r espectively of the V838 Mon outburst from |M unari| 

, .. ^ distance to 

extinction of 


VSljS Mon or b.l Koc 

Soarks et al.l 

2008k 

E(b-v) =0.7 (see Rus 

hton et al.|2U(J 

bk. 


that the peak magnitude was slightly brighter, although 
probably not significantly, given that the available pho¬ 
tometric data in all filters seems to suggest a smooth 
transition from the rise to decline phases of the light 
curve. 


3. SPECTROSCOPIC OBSERVATIONS 

We obtained four spectra of M31LRN 2015 using the 
SPRAT spectrograph on the LT in red optimized mode, 
with a resolution of 18 A, under photometri c conditions . 
A spectrophotometric standard (Hiltner 600; Stone 19771 
was observed on the night of 2015 February 4 in 1.4 see- 
ing (FWHM) and used to flux calibrate the data. SPRAT 
has a fixed width slit equivalent to 1'.'8 on the sky. See¬ 
ing in the science spectra varies in the range 1'.'5 — 2"5 
(FWHM) between nights. This seeing variation forms 
the major cause of uncertainty in our absolute flux cal¬ 
ibration, which we estimate as ^ 25%. However since 


all observations were obtained with the slit aligned with 
the parallactic angle, the relative flux calibration within 
a single spectrum is better constrained (typically < 10% 
between 5000 and 7500 A for all data). 

The first spectrum was obtained on 2015 January 
16.86, 3.2 days after discovery, and prior to peak. The 
only strong feature is the Ha emission, which has a 
FWHM of 900 ± 100 kms"! and peaks at 6551.8 ± 0.5 A, 
consistent with an M31 origin. We derive a Ha flux of 
(7.8 ±0.7) X 10“^^ ergcm“^ s“^ from the best-fit Gaus¬ 
sian profile. A small amount of Na i D absorption is de¬ 
tected (FWHM = 600±200kms“^), although this could 
be interstellar in origin, as well as H,5 being tentatively 
identified. This, along with the other SPRAT spectra 
taken of M31LRN 2015, is shown in Figure!^ 

The second spectrum was obtained 17 clays later on 
February 2.86. In this spectrum, the Ha flux had weak¬ 
ened to (3.6 ± 1.0) X 10“^^ ergcm“^ s“^, the continuum 
is much redder, Na i D absorption is now strong and 
several other strong absorption features have emerged. 
The third spectrum, taken February 12.84, is somewhat 
similar, but Ha is now only just visible above the contin¬ 
uum (which is even redder). The most prominent ab¬ 
sorption features in the second and third spectra are 
at 5882 ± 1 A (Na I D; FWHM = 1000 ± 200kms"!), 
6131 ± 3, 6240 ± 1, and 6482 ± 1 A. The 6482 ± 1 A ab¬ 
sorption line (FWHM = 900 ± 200kms“^) is likely to 
be Ba ii (6496.9 A) blueshifted by 690± 50kms“^, mak¬ 
ing the more poorly constrained 6131 ± 3 A line consis¬ 
tent with Ba II (6141.7A). There is also a less promi¬ 
nent absorption feature blue-ward of Na i D with a min¬ 
imum at 5841.6 A, but a profile could not be fit, which 
is possibly caused by Ba ii (5853.7 A). The absorption at 
6240 ± 1 A, may be caused by Fe ii, as seen ( along with 
Ba ii) in the V838 Mon spectra published by |Goranskii| 
(|2002|). We estimate the Na i D equivalent width 
as 1.2 ± 0.3 A on January 16.86 and 5 ± 1 A on 


et al. 

(EWd 


Febru ary 12.84. Using the relation from |Poznanski et al.| 
(20111, this gives a rough approximation of ~ 0.4 

on January 16.86 and E(^b-v) ~ 2 on February 12.84. 
Note however that those authors warn against using Na i 
D EW in low-resolution spectra to estimate absolute ex¬ 
tinction, but it does at least show that the absorbing 
column may be increasing between the two epochs. 

The fourth spectrum was taken February 28.85 and 
showed signihcant changes, with the spectrum now some¬ 
what resembling that of a K or M-type supergiant. No 
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Figure 2. Top: the four LT SPRAT spectra of M31LRN 2015 taken 2015 January 16.86, February 2.86, February 12.84 and February 
28.85. with peak occurring between 2015 January 19.89 and January 22.84 (see text). Bottom: three spectra of V838 Mon from |Rushton| 
|et al.| ||2005[l with respect to its peak, binned to a similar resolution as our SPRAT spectra of M31LRN 2015. 


significant Ha emission is detected and there are clear 
TiO absorption bands, with the Ba ii and Na i D lines 
still present. A spectrum was also taken on February 
25.84 but is not included as it was very similar to the 
one taken February 28, but noisier. 

4. THE PROGENITOR SYSTEM 


is undetected in the F555W data. The calculated posi¬ 
tion of M31LRN 2015, along with those nearby sources 
are shown in Figure Using a sim ilar method to that 

described in .. 

probability o: 


Williams et al. (2014), we determine the 
a source being within 0'1027 by chance is 


just 3%, given the resolved local stellar density. 


Having successfully r ecovered a number of progenitor 
systems for M31 novae (|B ode et al.||2009 Darnley et al. 
'2014| iWilliams et al.|2U14|), we used the same method to 
search for a possible progenitor system of M31LRN 2015. 
The archival Hubble Space Telescope (HST) observations 
were taken using the Advanced Camera for Surveys Wide 
Field Channel (ACS/WFC) through F555W and F814W 
filters on 2004 August 16.5 (prop. ID: 10273). Three 
F-band LT frames taken on 2015 January 17.87 were 
stacked and 37 bright stars common to both observations 
were used to compute the geometri c transform ation be¬ 
tween the two datasets using IRAF (Tody 19861 routines. 
This transformation was used to determine the position 
of the quiescent M31LRN 2015 in the archival HST data. 

This analysis reveals a source within 1 cr of the position 
of M31LRN 2015 in the HST data. Photometry was per¬ 
forme d on all sources using DOLPHOT (v2.0^; Dolphin 


2000). The closest source to the calculated position of 
A131LRN 2015, the progenitor candidate, is located at 
a distance of 0.55 ACS/WFC pixels, 0.5 (t or 0"027, and 
has magnitudes of F555W = 23.43 ± 0.04 and F814W 
= 22.09 ± 0.01 (these are in agreement with a similar 
analysis by Dong et al. 2015). We also note there is 
an additional source 1.4cr or 0'.'077 away. This second 
source has a magnitude of F814W = 24.48 ± 0.07, but 


^ http://purcell.as.arizona.edu/dolphot 



Figure 3. F814W HST (negative) image of the location of 

M31LRN 2015 taken 2004 August 16.5. The Ict and 3cr errors 
on the calculated position of the quiescent system are represented 
by green circles. The red crosses represent the positions of the 
nearby sources in the HST data. 


Emp loying the conversions from [Sirianni et al^ 
(2005) and using the distance and_reddening toward 
M31LRN 2015 in M31 (see Section we estimate the 
absolute magnitude of the progenitor candidate to be 
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Mv = -1.50 ± 0.23, Mi = -2.55 ± 0.13 with {V - I)o = 
1.05 ± 0.15, consistent with a red giant. 

The Local Grou p Galaxies Survey (LGGS; 
Massey et al. 2006) data reveal a nearby source. 


JUU42U8.U6+4055Orr at ^ = 22.33 ± 0.09, {B - V) 
= 0.75 ± 0.17, {V - R) = 0.62 ± 0.10 and (i? - I) 
= 0.65 ± 0.04. While it would first appear that the 
progenitor is brighter in thejMassey et al. (2006) catalog, 
an inspection of the images reveals this source is in 
fact a blend of several of the stars visible in Figure 
including the two brightest. A source is clearly visible at 
the position of M31LRN 2015 in the LGGS narrowband 
Ha data (taken 2002 September 11), but nothing is 
seen at this position in either the [O ill] or [S ii] data. 
To remove any contribution from continuum emission 
in these narrowband images we subtracted the scaled 
LGGS i?-band image from both the Ha and [S ii] 
images, and the scaled F-band image from the [O ill] 
image. The scaling was calculated by determining a 
linear fit between the broadband and narrowband data 


using SExtractor (v2.19.5; Bertin & Arnouts 1996) 
photometry of ^ 20,000 objects in each imag^ A 
clear source was seen in the continuum-subtracted Ha 
data, indicating a strong Ha exc ess in the spectrum o f 
the pre-outburst M31LRN 2015. Azimlu et al. (2011), 
who also used the LGGS data, list it as an H ii region 
(^1527 in their paper), with an estimated luminosity 
of 6 X lO^^ergs”^. If indeed this is a star forming 
H II region, the luminosity implies a star formation 
rate of ~ 5 x 10“^ Mq yr“^ (using the conversion from 


Kennicutt 1998). A SN remnant origin for this emission 
seems unlikely, due to the lack of [S ii] emission and 
the small size, which implies an age of less than a 
few hundred years. If this is not an H ii region, the 
most likely source of the Ha emission is the progenitor 
system itself, either directly from the progenitor or a 
companion, or from a period of mass loss. 

5. DISCUSSION 

M31LRN 2015 is a luminous, red transient in M31, 
characterized by weakening Ha emission on an increas¬ 
ingly red continuum, with a number of absorption lines 
emerging after peak, including Na i D and Ba ii, and 
TiO bands appearing at later times. Considering the 
spectroscopic and photometric evolution of the outburst 
we concl ude that this object is a LRN in M31, a conclu¬ 
sion that Kurtenkov et al. (2015b) reached independently 
from their data. We now compare the properties of this 
system to those of other proposed LRNe and ILRTs. 

The light curves of both LRNe and ILRTs are char¬ 
acterized by slow evolution (compared to most novae) 
and significant redde ning as they fade. The light curves 
of LRNe V1309 Sco ( [Mason et al.|[Mo| ) and V838 Mon 
(|Munari et al.|2002b|) show the brightness quickly falling 
in the bluer hlters, but plateauing in the red filt ers. The 
light curve of V4332 Sagittarii (another LRN; Martini 
et al.|199'9 ) also shows this object to be reddening about 
two weeks after discovery, although there is no color in¬ 
formation for the earlier periods. The only LRN that 
shows clear evidence of a significant multiple-peak light 
curve structure is V838 Mon. ILRTs are more luminous 
at peak than the LRNe discussed above and they also 
show different photometric evolution. The first ILRT 
discovered, M85 OT2006-1, shows a long plateau phase 


(^ 70 days), with relatively little color evolution un¬ 
til later times (|Kulkarni et al. 2007). ILRT RTF lOfqs 
showed a plateau phase (^ 40 days), before rapidly red¬ 


dening (Kasliwal et al. 
also s hows no ma.ior op 


2011), and NGG 300 OT2008 
ical color evolution until later 
times dBond et al.||2009|). Our observations show that 
M31LRiN 2015 clearly resembles LRNe, rather than IL¬ 
RTs, both in peak luminosity and light curve evolu¬ 
tion. In the late-time M31LRN 2015 light curve (see 
Figure [H, the F-band plateau appears to be ending, giv¬ 
ing a plateau time of 40—50 days. Comparing the plateau 
time and lumino sity to the common-e nvelope event theo¬ 
retical models of Ivanova et al. (|2013 see their Figure 1), 
M31LRN 2015 is consistent with the model for a merger 
with ejecta mass of 0.1 Mq, but is also consistent with 
an event where the entire envelope is ejected. 

The spectra of LRNe and ILRTs do show similarities, 
with both types of object exhibiting strong Na i D ab¬ 
sorption and narrow Ha emission, but they also appear 
to show some differences. Strong [Ca ii] emission ha s 
been seen in several ILRTs (see e.g. Kasliwal et al.|2011 ), 
which is not the case for LRNe. Additionally, strong 'I'iO 
abso rption bands were s een at later times in LR Ne V4332 
Sgr dMarti ni et al.|1999|) and V838 Mon (see e. g.lRushton 
et al.|2UU5|). I'hese were also seen in M31-RV (Rich et al. 


1989[ ). V838 Mon is the LRN with the best spectroscopic 


coverage, and the evolution of M31LRN 2015 seems to 
match that very well, for example, Ba ii (6496.9 A) ap¬ 
pears between ~ 2 days prior to peak (t ~ —2 days; on 
2002 February 5) and t ^ +9 days in V838 Mon (|Goran- 
skii et al.|2002 ), and also appears in M31LRN 2015 after 


peak. It can be seen from Figure]^ that our second and 
third spectra of M31LRN 2015 are remarkably similar to 
spectra t aken around t ~_+32 days of V838 Mon pub¬ 
lished by|Rushton et al.|(|2005|), with our fourth spectrum 
of M31LRiN 2015 similar to V838 Mon at t ^ -1-60 days. 

The only significant differences between M31LRN 2015 
and V838 Mon are the pre-outburst detections of the pro¬ 
genitor systems and the second bright peak in the light 
curye of the latter. The V838 Mon syst em contains a 
young main-sequence binary companion (Munari et al. 


2002a). For the case of M31 LRN 2015, we may again be 


seeing a binary companion, but at a more eyolyed stage 
and depending on the mechanism, this may not repre¬ 
sent different progenitor stars for the outbursts them- 
selyes. Notably, no young stellar population similar to 
the locale of V838 Mon was fo und to be coincident with 
M31-RV ( Bond &: Siegel|[2006 ). If the giant planet cap¬ 
ture mechanism was responsible for LRNe outbursts, this 
system (along with M31-RV) would represent the first 
(albeit indirect) detection of extragalactic planetary sys¬ 
tems. While the detected progenito r candidate does not 
exclude a classical noya mechanism (Darnley et al.|[^12 


Williams et al.|[2014 ) for M31LRN 2U15, the presence o: 


a high-mass mam-sequence star in the V838 Mon system 
and the strong eyidence for a stellar merger in the case 
of V1309 Sco appears to rule it out for those systems. 
The spectra ar e not consistent wit h a typical noya erup¬ 
tion, although Shara et al. (2010) offer an explanation 
for this. The strong eyidence tor the V1309 Sco outburst 
being cause d by the merger of a contact binary (Tylenda 


et al. 2011) suggests that if a single progenitor mecha- 


nism IS responsible for all LRNe, this is li kely to be the 
stellar merger scenario. Dong et al. (2015) reported pos- 
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sible pre-outburst variability in M31LRN 2015 from dif¬ 
ference image analysis (DIA) of archival data. Although 
the presence of many faint sources in such crowded fields 
(see Section can affect the reliability of photometry, 
even when employing DIA, it is unlikely in this case and 
the significant brightening they detect from g ~ 22.6 to 
g ^ 20.8 between 2009 and 2014 certainly warrants fur¬ 
ther investigation. 

6. CONCLUSIONS 

In this Letter we have presented photometric and spec¬ 
troscopic observations of M31LRN 2015, and a search for 
the progenitor system in broad and narrowband archival 
data. Here we summarize the main conclusions. 

1. We used the LT to photometrically and spectro¬ 
scopically observe M31LRN 2015, which is a LRN 
peaking at an absolute magnitude of My = —9.4± 
0.2. It took several days to reach optical peak, 
with the optical colors becoming increasingly red¬ 
der over the following weeks. 

2. The spectra are characterized by a reddening con¬ 
tinuum, weakening Ha emission and emerging ab¬ 
sorption features, including Na i D and Ba ii. At 
later times TiO absorption bands are also seen. 
M31LRN 2015 shows a number of similarities with 
the spectra of other transients classified as LRNe, 
resembling V838 Mon particularly closely. 

3. The likely progenitor system has an absolute mag¬ 
nitude of My = —1.50 ± 0.23 with {V — I)q = 
1.05 ± 0.15, although this could represent photom¬ 
etry of a more luminous companion star. Pre¬ 
outburst Ha emission is likely to be related to the 
progenitor system, either directly or the local envi¬ 
ronment. 

4. The nature of all of the LRN class of objects is still 
uncertain, but our observations of M31LRN 2015 
and its likely progenitor system suggest that the 
mechanism(s) must produce very similar outbursts 
(in peak luminosity, timescale and outburst evolu¬ 
tion), and possibly within differing stellar popula¬ 
tions. If a single progenitor mechanism is respon¬ 
sible for all LRNe, the balance of evidence favors 
the stellar merger pathway. 

Further observations, particularly of the post-outburst 
system, may yield further clues to the nature of 
M31LRN 2015 and the LRN class of objects in general. 

We would like to thank M. Rushton for supplying 
the V838 Mon spectra used here for comparison with 
M31LRN 2015. We also thank an anonymous referee 
for the constructive suggestions. The Liverpool Tele¬ 
scope is operated on the island of La Palma by Liverpool 
John Moores University in the Spanish Observatorio del 
Roque de los Muchachos of the Institute de Astrofisica 
de Canarias with financial support from the UK Science 
and Technology Facilities Council. Some of the data 
presented in this paper were obtained from the Multi¬ 
mission Archive at the Space Telescope Science Institute 


(MAST). STScI is operated by the Association of Uni¬ 
versities for Research in Astronomy, Inc., under NASA 
contract NAS5-26555. Support for MAST for non-HST 
data is provided by the NASA Office of Space Science via 
grant NNX09AF08G and by other grants and contracts. 
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